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ABSTRACT 

We analyze the small-scale clustering in "MegaZ-LRG" , a large photometric-redshift 
catalogue of Luminous Red Galaxies extracted from the imaging dataset of the Sloan 
Digital Sky Survey. MegaZ-LRG, presented in a companion paper, spans the redshift 
range 0.4 < z < 0.7 with an r.m.s. redshift error ~ 0.03(1 -I- z), covering 5,914 
deg^ to map out a total cosmic volume 2.5 Gpc'^. In this study we use 380,000 
photometric redshifts to measure significant deviations from the canonical power-law 
fit to the angular correlation function in a series of narrow redshift slices, in which we 
construct volume-limited samples. These deviations are direct signatures of the man- 
ner in which these galaxies populate the underlying network of dark matter haloes. We 
cleanly delineate the separate contributions of the "1-halo" and "2-halo" clustering 
terms and fit our measurements by parameterizing the halo occupation distribution 
N{M) of the galaxies. Our results are successfully fit by a "central" galaxy contribu- 
tion with a "soft" transition from zero to one galaxies, combined with a power-law 
"satellite" galaxy component, the slope of which is a strong function of galaxy lumi- 
nosity. The large majority of galaxies are classified as central objects of their host dark 
matter haloes rather than satellites in more massive systems. The effective halo mass 
of MegaZ-LRG galaxies lies in the range logj^QiM^g/h-^MQ) = 13.61 13.80 (in- 
creasing with redshift, assuming large-scale normalization erg — 0.8) for corresponding 
number densities in the range Ug = 5.03 0.56 x lO^"' Mpc^^. Our results confirm 
the usefulness of the halo model for gaining physical insight into the patterns of galaxy 
clustering. 
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1 INTRODUCTION 

Photometric redshift surveys offer a route to delineating 
the large-scale structure of the Universe that is increasingly 
competitive with spectroscopic redshift surveys (Budavari 
et al. 2003; Seo & Eisenstein 2003; Amendola, Quercellini 
& Giallongo 2004; Dolney, Jain & Takada 2004; Blake & 
Bridle 2005; Phleps et al. 2006; Zhan et al. 2006, Blake et 
al. 2007; Padmanabhan et al. 2007). The ease with which 
modern imaging surveys can map large areas of sky to faint 
magnitude limits compensates for the absence of precise 
(but time-consuming) spectroscopic redshift measurements 
for individual galaxies. An absolute pre-requisite, however, is 
the availability of high-quality photometric galaxy redshifts 
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with known error distributions (established for example via 
spectroscopy of sub-samples), together with accurate survey 
photometric calibration over large angles of sky. Recent ob- 
servational efforts have enabled both of these criteria to be 
satisfied. 

The Sloan Digital Sky Survey (SDSS; York et al. 2000) 
has now provided an accurately-calibrated imaging dataset 
over roughly a fifth of the sky, which can be used to extract 
samples of galaxies in a uniform manner. In particular, a 
photometric catalogue of Luminous Red Galaxies (LRGs) 
can be readily extracted using a series of well-understood 
colour and magnitude cuts (Eisenstein et al. 2001). Owing to 
their high luminosity and typical residence in the most mas- 
sive dark matter haloes, LRGs are efficient tracers of cosmic 
structure across large volumes (Brown et al. 2003; Zehavi et 
al. 2005a; Eisenstein et al. 2005b). Moreover, these objects 
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provide particularly reliable photometric redshifts owing to 
the strong spectral break at « 4000A in the galaxy rest 
frame and the consequent rapid variation with redshift of 
their observed colours in the SDSS filter system (Padman- 
ablian ct al. 2005; CoUistcr ot al. 2007). Furthermore, the 
photometric redshift error distribution of the LRGs can be 
accurately calibrated owing to the existence of spectroscopic 
observations of a sub-sample as part of the 2dF-SDSS LRG 
and Quasar (2SLAQ) survey at the Anglo-Australian Tele- 
scope (Cannon et al. 2006). 

The combination of these datasets has allowed us to 
construct a large catalogue of more than 10® LRGs spanning 
the redshift range 0.4 < z < 0.7 with an r.m .s. photometric 
redshift error = a/< (Sz)^ > - < Sz >^ ^ 0.03(1 -|- z) 
where Sz = Zphot — Zspec- We have dubbed this catalogue 
"MegaZ-LRG" (Colhster et al. 2007). The database covers 
almost 6,000 deg^, or an effective volume « 2.5 Gpc^. 
Wo have already used this catalogue to measure the large- 
scale clustering of the galaxies on linear and quasi-linear 
scales via a power spectrum analysis and thereby extract 
measurements of cosmological parameters (Blake et al. 2007; 
see also Padmanabhan et al. 2007). In this study we turn to 
the small-scale clustering properties of the LRGs. Our goal is 
to present new measurements of the small-scale correlation 
function of LRGs at 2 ~ 0.5, and to connect these measure- 
ments to the physical manner in which these LRGs populate 
dark matter haloes. 

It has long been known that different classes of galaxy 
possess different clustering properties, in a manner con- 
nected to their small-scale environments (for recent observa- 
tional studies we refer the reader to Norbcrg et al. 2002; Bu- 
davari et al. 2003; Hogg et al. 2003; Zehavi ct al. 2005b; and 
references therein). For many years these differing clustering 
properties could be adequately described by fitting a simple 
power-law function to the two-point galaxy correlation func- 
tion on small scales (Peebles 1980). However, recent surveys 
have measured the clustering pattern accurately enough to 
detect deviations from the canonical clustering power-law 
(e.g. Hawkins et al. 2003; Zehavi et al. 2004; Zheng 2004; 
Eisenstein et al. 2005a; Zehavi et al. 2005a,b; Phleps et al. 
2006). These deviations provide an important insight into 
the processes of galaxy formation. 

The richer structure in the galaxy clustering pattern re- 
vealed by recent surveys has been successfully interpreted 
in terms of the "halo model" (e.g. Seljak 2000; Peacock 
& Smith 2000; Scoccimarro et al. 2001; Cooray & Sheth 
2002; Berlind & Weinberg 2002; Kravtsov et al. 2004; Ze- 
havi et al. 2004; Zheng 2004; Zheng et al. 2005; Zehavi et al. 
2005b; Colhster & Lahav 2005; Tinker et al. 2005). In this 
model, the small-scale clustering of a distribution of galaxies 
is linked to the underlying network of dark matter haloes, 
whose properties can be measured using cosmological sim- 
ulations. A class of galaxies is assumed to populate haloes 
in accordance with a statistical "halo occupation distribu- 
tion" as a function of the halo mass. The clustering then 
naturally separates into two components: the distribution of 
galaxies within individual haloes, which dominates on small 
scales ( < 1 Mpc) , and the mutual clustering of galaxies in- 
habiting separate haloes, which dominates on larger scales 
( > 1 Mpc). The combination of these two terms can accu- 
rately model the observed scale-dependent features in the 
small-scale clustering pattern. 



The outline of this paper is as follows: in Section 2 we 
briefly describe the "MegaZ-LRG" dataset. We then present 
measurements of the angular correlation function in narrow 
redshift slices together with simple power-law fits in Section 
3. In Section 4 we introduce the framework of the halo model 
and in Section 5 we re-fit the clustering measurements by pa- 
rameterizing the halo occupation distribution of the LRGs, 
comparing our results to previous work in Section 6. Section 
7 investigates a range of potential systematic photometric 
errors in the catalogue that may bias our results. We con- 
clude in Section 8. 

Throughout our study we assume a fixed set of large- 
scale cosmological parameters: fractional matter density 
Qm = 0.3, cosmological constant Qa = 0.7, curvature 
fik = 0, Hubble parameter h = 0.7, fractional baryon den- 
sity /b = fib/f^m = 0.15, slope of the primordial power 
spectrum rig = 1, and overall normalization of the power 
spectrum as — 0.8. These values are consistent with fits 
to the large-scale power spectrum of the LRGs (Blake et 
al. 2007) and to the latest measurements of the anisotropy 
spectrum of the Cosmic Microwave Background (Spergel et 
al. 2007). We note that there is some degeneracy between 
the assumed cosmological parameters and the fitted halo 
model parameters, and we investigate how the best-fitting 
halo model parameters depend on the values of as and Ub- 



2 THE DATA SET 

We analyze angular clustering in the "MegaZ-LRG" galaxy 
database, a photometric-redshift catalogue of Luminous Red 
Galaxies (LRGs) based on the imaging dataset of the SDSS 
4th Data Release. The construction of this catalogue is de- 
scribed in detail by Colhster et al. (2007), and cosmological- 
parameter fitting to the angular power spectrum is presented 
by Blake ct al. (2007) . We only provide a brief description of 
the catalogue here, refering the reader to these two papers 
for more information. 

MegaZ-LRG contains over 10® LRGs spanning the red- 
shift range 0.4 < z < 0.7 with an r.m.s. redshift error 
5z ~ 0.03(1 + z). The angular selection function is de- 
scribed by Blake et al. (2007) and encompasses 5914 deg^ 
(the three southern SDSS stripes are excluded). The sam- 
ple was selected from the SDSS imaging database using a 
series of colour and magnitude cuts (Eisenstein et al. 2001; 
Colhster et al. 2007), which amount to a magnitude-limited 
sample of Luminous Red Galaxies with i-band magnitudes 
17.5 < i < 20. Photometric redshifts were derived for 
these galaxies using an Artificial Neural Network method, 
ANNz (Firth, Lahav & Somerville 2003; Colhster & La- 
hav 2004), constrained by a spectroscopic sub-sample of 
« 13,000 galaxies obtained by the 2SLAQ survey (Cannon 
et al. 2006). In this paper we analyze a "conservative" ver- 
sion of the database in which the selection cuts applied to 
the imaging database are identical to those used to produce 
the great majority of the spectroscopic sub-sample (e.g. the 
faint magnitude limit is brightened to i = 19.8). Star-galaxy 
separation cuts were applied both in the initial selection 
from the SDSS database and via the neural network anal- 
ysis. These cuts produced a catalogue with 644,903 entries, 
amongst which there is a 1.5% M-star contamination (see 
Blake et al. 2007). 
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We used the photometric redshifts to divide the sam- 
ple into four narrow redshift slices of width A^phot = 0.05 
between z = 0.45 and z = 0.65 (there are very few galax- 
ies in the catalogue outside this redshift range). We then 
applied a luminosity threshold in each redshift range to 
create a "volume-limited" sample of gala^Kies, as assumed 
by the halo model. The luminosity threshold is given by 
the faint apparent magnitude limit i = 19.8 applied at 
the most distant redshift of each slice. We calculated lu- 
minosities for each galaxy using Luminous Red Galaxy K- 
corrections from the 2SLAQ survey (Wake et al. 2006), 
not including an evolution correction. The faint absolute 
i-band magnitude limits for each redshift slice arc then 
Mi-51ogi(|/i = (-22.23, -22.56, -22.87, -23.20). The rmm- 
ber of galaxies remaining in each redshift slice after the lu- 
minosity cut was N = (168287, 118863, 70229, 27203) with 
corresponding surface densities (28.5,20.1,11.9,4.6) deg~^. 
The spectroscopic redshift distribution of galaxies in each 
photo-z slice (including the luminosity cut) can be deduced 
using the 2SLAQ spectroscopic sub-sample. As shown by 
Blake et al. (2007), the spec-z probability distribution for 
each slice is well-described by a Gaussian function; the mean 
/i and standard deviation a for each photo-z slice are listed 
in Table 1 (note that these values are slightly different from 
those listed in Blake et al. 2007, due to the additional lumi- 
nosity threshold applied in the current study). These red- 
shift distributions are used to project the model correlation 
function to fit the observed angular clustering in each slice. 



3 ANGULAR CORRELATION FUNCTION 
MEASUREMENTS 

3.1 Method 

We used the Landy & Szalay (1993) estimator to mea- 
sure w{9) for each photometric redshift slice in 30 
logarithmically-spaced angular separation bins between 6 = 
0.001° and 9 = 1°. For each redshift slice we generated 10 
random datascts across the survey geometry, each contain- 
ing the same number of galaxies as the survey datasets. We 
assume that the 1.5% stellar contamination is distributed 
evenly across the redshift slices in an unclustered fashion 
such that the amplitude of the measured angular correla- 
tion function is simply reduced by a constant factor (1 — f)'^ 
where / = 0.015. We corrected our estimates of w{9) and 
corresponding errors upwards by this factor (of 3.0%). 



3.2 Error determination 

We estimated the covariance matrix of the errors in the sepa- 
ration bins using the technique of jack-knife resampling. For 
each measurement of w{9) we divided the survey area into 
N = 393 sub-fields of constant area (of 13.3 deg'^) using a 
grid of right ascension and declination divisions (see Figure 
1). When creating the grid, we allowed a sub-field to contain 
up to 20% fractional area beyond survey edges or of survey 
holes. The rmmbor of sub-fields was chosen to be sufficient 
for estimating each unique element of the covariance ma- 
trix used in the model fitting with statistical independence. 
However, we checked that our best-fitting parameters did 
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Figure 1. The grid oi N = 393 sub-fields of equal area, defined by 
constant right ascension and declination boundaries, which was 
used to estimate the error in the angular correlation function via 
jack-knife re-sampling of the measurement. 



not depend on the number of sub-fields, or on the restric- 
tion of our fits to the most significant principal components 
of the covariance matrix using singular value decomposition 
techniques. 

Having defined the grid of sub-fields, wc then mea- 
sured w{0) N times, which wc label as Wi{9) from i = 1 
to i = N, in each case omitting just one of the sub-fields 
(and using the remaining N — 1 fields). The covariance of 
the measurements between separation bins j and k was de- 
duced as: 



Cjk = {w{9j) wiOk)) - {w{ej)) {w{9k)) 



(1) 



(iV- 1) I ^ W{9j) X w{9k) 1(2) 



where w{9,) =J2^ wA9,)/N. 

Figure 2 plots the jack-knife errors in the separation 
bins for each redshift slice, normalized by the error deter- 
mined assuming simple Poisson statistics (for which the er- 
ror in the data pair count of DD objects in a separation 
bin is V DD). For most of the angular range under investi- 
gation, the jack-knife errors are less than 50% higher than 
those predicted by Poisson statistics. For the largest angular 
scales 9 considered, the jack-knife error increases relative to 
the Poisson error, owing to the increasing importance of edge 
effects and the heightened "cosmic variance" owing to the 
reduced number of independent cells of size 9 that can be ac- 
commodated by the dataset. This is the familiar result from 
clustering measurements that Poisson noise dominates on 
small scales, and cosmic variance dominates on large scales. 
The full covariance matrices are displayed in Figure 3 by 
plotting in grey-scale the "correlation coefficient" between 
two separation bins i and j: 



Cov(i, j) 



^yCov{^,i) Cov{j,j) 



(3) 



Angular correlation function measurements in largo separa- 
tion bins are positively correlated. The amplitude of these 
correlations decreases with redshift because, as the number 
density of the sample reduces, Poisson noise becomes more 
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Figure 2. The error in the angular correlation function deter- 
mined by jack-knife re-sampling for the four redshift slices, nor- 
malized by the result assuming simple Poisson statistics. 
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Figure 3. Grey-scale plot of the correlation coefficient r of equa- 
tion 3, indicating the degree of covariance between different sepa- 
ration bins for each redshift slice. The correspondence of r to the 
grey-scale is chosen as r = —0.5 (white) to r = 0.5 (black). 



importamt compared to cosmic variance. These covariance 
matrices are always used in our model fitting. 



3.3 Power-law fits 

As an initial step we fitted power-laws w{0) = aO^ '' to 
the measured angular correlation functions in the four red- 
shift slices (using clii-squared minimization with the full co- 
variance matrix). We excluded the first 3 data points with 
separations d < 0.002° (= 7 arcsec), as these scales are po- 
tentially afi^ected by astronomical seeing and issues of galaxy 
merging and blending (the decrease in the value of the cor- 
relation functions at the smallest scales in Figure 4 is not 
a property of galaxy clustering). The best-fitting power- law 
parameters (0,7) are listed in Table 1 and the models and 



data points are displayed in Figure 4. The best-fittiug slopes 
7 for the redshift slices take values in the range 1.94 1.96, 
consistent with previous studies of Luminous Red Galaxies 
(e.g. Eisenstein et al. 2005a; Zehavi et al. 2005a). However, 
as indicated by the high values of the minimum clii-squared 
statistic Xpi in Table 1 compared to the 25 degrees of free- 
dom, a power-law is not a good fit to the data. The fact that 
we detect deviations from simple power-law clustering with 
high significance indicates that our photo-2 survey cam be 
used to fit more complex and physically-insightful models to 
the small-scale clustering pattern. 

Nevertheless, we can use the power-law amplitude of the 
angular clustering to estimate the three-dimensional cluster- 
ing length ro of the population of galaxies, using the known 
redshift distribution p{z) in each slice. If wo define the spa- 
tial two-point correlation function ^(r) = (r/ro)"'' then 
the amplitude of the power-law angular correlation function 
w{0) = aO^ '' follows from Limber's equation: 



a 



dzp{i 



2 /cte 
' \dz 



x{zf~ 



(4) 



where x[z) is the co-moving radial co-ordinate at redshift 
z and Ct- = r(i)r(2 - i)/r(^). In Table 1 we determine 
the corresponding values of ro and 7 for each redshift slice. 
The amplitude of the clustering length, ro = 7.2 — » 9.3 
h^^ Mpc, is consistent with highly-biased massive galaxies. 
Our small-scale angular clustering measurements are not af- 
fected by redshift-spaco distortions, since the photo- 2: errors 
are much bigger than the peculiar velocities of the galax- 
ies. Correlated redshift-space distortions are important on 
larger scales however, as discussed in Blake et al. (2007) 
and Padmanabhan et al. (2007). 

The clustering amplitude systematically increases with 
redshift for two reasons: 

(i) Each volume-limited sample in successive redshift 
slices has the same limiting apparent magnitude hence 
higher luminosity threshold. The highor-rodshift galaxies are 
hence preferentially more luminous and more strongly clus- 
tered (e.g. Norberg et al. 2002; Zehavi et al. 2005b). In Ta- 
ble 1 we list the threshold absolute i-band magnitude Mi of 
galaxies in each redshift slice, calculated using a Luminous 
Red Galaxy K-correction (Wake et al. 2006). 

(ii) In standard models of the evolution of galaxy clus- 
tering, the bias factor of a class of galaxies increases with 
redshift in opposition to the decreasing linear growth fac- 
tor, in order to reproduce the observed approximate con- 
stancy of the small-scale co-moving clustering length (e.g. 
Lahav et al. 2002). Simple models for this effect such as 
h{z) = l + {ho- l)/D{z) (Fry 1996) or h{z) = bo/D{z), 
where D{z) is the linear growth factor, predict an evolution 
in bias across our analyzed redshift range of A6 w 0.2. 

These trends are in good agreement with our measurements 
of the amplitude of the large-scale clustering pattern (Blake 
et al. 2007). 



4 HALO MODEL FRAMEWORK 

We used the halo model of galaxy clustering to produce 
model spatial correlation functions ^(r) to fit to our mear- 
surements. We summarize the ingredients of our model here. 
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Figure 4. The measured angular correlation function in the four redshift slices together with the best-fitting power laws. The turnover 

of the correlation function towards very small angular separations is not a property of galaxy clustering; wc restrict our fit to separations 
6 > 0.002°. The co-moving spatial scale corresponding to angular separation 0.1° at the median redshift of each slice is displayed in the 
bottom left-hand corner of each plot. 



Table 1. Power-law fits to the angular correlation function in the four redshift slices. Column 2 records the surface density of the galaxy 
sample in each redshift slice. Column 3 is the threshold galaxy absolute i-band magnitude in each slice for the volume-limited sample. 
Columns 4 and 5 list the parameters of the Gaussian redshift distribution for each slice, p{z) oc exp [—{z — /u)^/2<t^]. Columns 6 and 7 
show the best-fitting power-law parameters w{9) = aO^ where 6 is in degrees. The associated errors for a and 7 were obtained by 
marginalizing over values of the other parameter. Note that the amplitudes have been corrected upwards by a stellar contamination factor 
(1 — f)~^ = 1.03. Column 8 is the corresponding chi-squarcd statistic Xpi evaluated using the full covariancc matrix (which wc compare 
with 25 degrees of freedom). Columns 9 lists the inferred galaxy clustering length ro of the spatial correlation function ^(r) = (r/ro) ' ■ 



Redshift slice 


Density 


Mi -5 logio/i 




cr 


10^ X a 


7 


Xpi 


ro 




(cleg-2) 








(dcg^-i) 




(/i-i Mpc) 


0.45 < 2 < 0.5 


28.5 


-22.23 


0.476 


0.034 


2.84 ±0.12 


1.96 ±0.01 


177.1 


7.2 ±0.2 


0.5 < z < 0.55 


20.1 


-22.56 


0.528 


0.040 


2.79 ±0.14 


1.96 ±0.01 


132.8 


8.0 ±0.2 


0.55 < 2 < 0.6 


11.9 


-22.87 


0.574 


0.043 


2.93 ±0.18 


1.94 ± 0.02 


71.8 


8.5 ±0.3 


0.6 < z < 0.65 


4.6 


-23.20 


0.629 


0.052 


2.78 ±0.19 


1.94 ± 0.02 


34.9 


9.3 ±0.3 



Further details can be found in e.g. Seljak 2000; Cooray & 
Sheth 2002; Berlind & Weinberg 2002; Kravtsov et al. 2004; 
Zchavi ct al. 2004; Zheng 2004; Zehavi et al. 2005b; Tinker 
et al. 2005. 

In the halo model framework, the clustering functions 

are expressed as a sum of components due to pairs of galaxies 
within a single dark matter halo (the "1-halo term" ^1) and 
to pairs of galaxies inhabiting separate haloes (the "2-halo 
term" ^2): 

^(r) = [l+6(r)l + 6(r) (5) 



where the "1+" at the start of the expression arises because 
the total number of galaxy pairs (oc 1 + C) is the sum of 
the number of pairs from single haloes (oc 1 + ^i) and from 
different haloes (oc 1 + ^2). The two terms dominate on dif- 
ferent scales, with the 1-halo term only important on small 
scales < 1 Mpc. 

The fundamental ingredient of the galaxy halo model 
is the halo occupation distribution (HOD), which describes 
the probability distribution for the number of galaxies N 
hosted by a dark matter halo as a function of its mass M. In 
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order to construct the 1-halo and 2-halo two-point clustering 
terms, we require the first and second factorial moments of 
the HOD, < N\M > and < N{N - 1)\M >. We make 
the assumption that the first galaxy to be hosted by a halo 
lies at the centre of the halo, and any remaining galaxies 
are classified as "satellites" and distributed in proportion 
to the halo mass profile. We apply different HODs for the 
central and satellite galaxies, < Nc\M > and < 7Vs|M > 
respectively, where 



< N\M >=< Nc\M > (1+ < N,\M >) 



(6) 



Equation 6 takes its form because a halo can only host a 
satellite galaxy if it already contains a central galaxy. We will 
use the notation Nc{M) =< Nc\M >, Ns{M) =< Ns\M > 
and N{M) =< N\M > in the equations that follow. 



4.1 The 1-halo term ^i(r) 

The 1-halo galaxy correlation function is composed of con- 
tributions from central-satellite pairs and satellite-satellite 
pairs. It is convenient to evaluate these two contribu- 
tions separately. The 1-halo correlation function for central- 
satellite pairs is given by: 



1 + 6 



I 



dMn{M) 



Nc{M)N,(M) p(r|M) 



ni/2 



M 



(7) 



where rig is the galaxy number density, n(M) is the halo 
mass function, and p{r\M) is the halo density profile. The 
lower limit for the integral is the halo mass M correspond- 
ing to a virial radius r, given that less massive haloes have 
smaller radii and cannot contribute any central-satellite 
galaxy pairs with co-moving separation r: 



(8) 



where p = 2.78 x 10"f^m ft^M© Mpc~^ is the co-moving 
background density of the Universe, and A = 200 is the 
critical overdensity for virialization. 

It is simplest to evaluate the 1-halo correlation function 
for satellite-satellite pairs in Fourier space (where convolu- 
tions become multiplications). The power spectrum is: 

= rdMniM)^^^iM}^Hk\M)f (9) 



where u{k\M) is the Fourier transform of the halo den- 
sity profile p{r\M). Because satellite galaxies arc Poisson- 
distributed, we can write < N^iN^ - 1) >=< A^^ >^ to 
obtain the above equation. The correlation function corre- 
sponding to equation 9 is then 



1 

2^ 



dkPi. 



sinfer 
kr 



(10) 



The total 1-halo correlation function is then derived as 

6 =6,c-s + Ci,.-<. (11) 



P2{k,r): 

P2{k,r) = Pm{k)x 



I 

JO 



dM n(M) b(M, r) ^^^} u{k\M) 
ng{r) 



(12) 



where Pm(k) is the non- linear matter power spectrum at 
the survey redshift, 6(M, r) is the scale-dependent halo bias 
at separation r, and n'g{r) is the restricted galaxy number 
density at separation r, where 

r-M„„(r) 



I'gir) 



dMn{M) N{M) 



(13) 



The mass truncation Miini('") must be included to incor- 
porate the effects of halo exclusion: more massive haloes 
would overlap at separation r. We derive the mass limit us- 
ing the "rig-matched" approximation described in Tinker et 
al. (2005). Firstly we calculate the restricted number den- 
sity using equation B13 in Tinker et al., which includes the 
effects of triaxiality: 
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dMin{Mi)N{Mi)x 



dMi n(M2) N{M2) P{r, Mi, Ma) 



(14) 



where P{r, Mi, AI2) quantifies the probability of non- 
overlapping haloes of masses Mi and M2 with separation 
r. Defining x = r/{Ri + R2), where Ri and R2 are the 
virial radii corresponding to masses Mi and M2, and using 
y = (a;-0.8)/0.29, then Tinker et al. obtain P{y) = 3y^-2y^ 
from simulations. Given this value of n'g (r) , we increase the 
value of Miini(r) in equation 13 to produce a matching num- 
ber density. This value is then used in equation 12 to produce 
the 2-halo power spectrum P2{k,r). 

Following Tinker et al. (2005), we assumed the following 
model for the scale-dependent bias: 



6^(M,r) = 6^(M) 



[l + 1.17g„(r)]^-- 



(15) 



[H-0.69 6n(r)]'-°' 

where (,m{r) is the non- linear matter correlation function. 
We derived matter power spectra using the "CAMB" software 
package (Lewis, Challinor & Lasenby 2000), including cor- 
rections for non-linear growth of structure using the fitting 
formulae of Smith et al. (2003) ("halofit=l" in CAMB). We 
outputted power spectra at the mean redshift of each slice, 
and obtained the correlation function using a Fourier trans- 
form. The 2-halo galaxy correlation function is obtained via 
the Fourier transform of equation 12: 



^2ir) 



•2^' Jo 



dkP2{k,r) k- 



sin kr 
kr 



(16) 



with the number of galaxy pairs corrected from the re- 
stricted galaxy density to the entire galaxy population: 

2 



1 + 6(0 = 



Tin 



[l + 52(r)] 



(17) 



4.2 The 2-halo term ^2(r) 

The 2-halo galaxy correlation function at separation r is 
evaluated from the scale-dependent 2-halo power spectrum 



4.3 Halo mass and bias functions 

The halo mass function n(M) describes the number density 
of haloes as a function of mass M. We introduce the new 
mass variable v following Press & Schechter (1974): 
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n{M)dM=j^f{v)dv 



(18) 



The new mass variable is defined hy p = [Ssc/i^iM, z)]^ , 
where Ssc is the hnear-thcory prediction for the present-day 
overdeusity of a region which would undergo spherical col- 
lapse at redshift z, and cr^{M, z) is the variance of the linear 
power spectrum in a spherical top hat which contains aver- 
age mass M: 



(19) 



where W{x) = (3/x^)[sina: — a: cos a:], M = ^nR^p, D{z) is 
the linear growth factor at redshift z, and f\in(A;) is the 
linear power spectrum at redshift zero. We approximate 
5sc = 1.686 independently of redshift. 

We use the Jenkins et al. (2001) model for the mass 
function: 



-oiexp[-|ln(^/5sc) +a2r^l 



(20) 



where ai = 0.315, a2 = 0.61 and 03 = 3.8. 

The halo bias function b{M) describes the biasing of 
a halo of mass Al with respect to the overall dark matter 
distribution. We use the Sheth, Mo & Tormen (2001) model 
for the bias function, with the revised parameters stated in 
Tinker at al. (2005): 



6(^) 



:l+;j-X 

qv + s{qv)^~* - 



-1/2 



(21) 



0.8. 



l + .(l-i)(l-|)(9z.)-* 
where the constants q = 0.707, s = 0.35 and t = 

4.4 Halo profiles 

We use the Navarro, Prenk & White (1997) dark matter halo 
density profile: 



p(r) = 



Ps 



(r/r,)(l-hr/r,)2 



(r < Tvir) 



(22) 



where is the chaxacteristic scale radius and pa provides the 
normalization. The profile is truncated at the virial radius 
Tvir, which is obtained from the halo mass via 



- f \ 



1/3 



(23) 



We parameterize the profile in terms of the concentration 
parameter c = rvir/rs- The normalization for the mass M is 



M 



p{r)Airr^dr = Anparl 



1 + c 



(24) 



We assume that the concentration parameter c depends on 
halo mass M and redshift 2; in a manner calibrated by nu- 
merical simulations (Bullock et al. 2001; Zehavi et al. 2004): 



c{M,z) 
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z \MqJ 



(25) 



where Mo is obtained from equation 19 by setting 
cr(Mo,0) = Ssc- For our adopted cosmological model we ob- 
tain Mo = 12.64 /i~^M0, resulting in a concentration c = 5 



for a halo of mass M = lO^^'/i^^M© at a = 0.5. This as- 
sumption is consistent with the measurement by Mandcl- 
baum et al. (2006b) of the concentration parameter for LRGs 
using galaxy-galaxy weak lensing. Fitting the concentration 
c as an extra free parameter produced no improvement in 
the minimum value of . In fact, there is a significant de- 
generacy between c and (3. 



4.5 Halo occupation distribution 

The halo occupation distribution (HOD) is a parameterized 
description of how galaxies populate dark matter haloes as 
a function of the halo mass M. As described above, we sepa- 
rate the distribution into separate HODs for central galaxies 
and for satellite galEixies. We adopt simple models motivated 
by results from simulations and semi-analytic calculations 
(e.g. Kauffmann, Nusscr & Steinmetz 1997; Benson et al. 
2000; Bcrlind et al. 2003; Kravtsov et al. 2004). 

For central galaxies, our basic model is a step function 
such that haloes above a minimum mass threshold Mcut con- 
tain a single central galaxy and haloes below this threshold 
contain no galaxies. However, our observational data con- 
sists of galaxy luminosities (with a scatter resulting from 
the photo-2 errors) rather than galaxy masses. We therefore 
follow Zheng et al. (2005) and "soften" the transition from 
zero to one galaxies with a further parameter Ucut such that 



< Nc\M >= 0.5 



1-l-erf 



logio(M/Mc, 

<7cut 



(26) 



The scatter in halo mass for a fixed galaxy luminosity results 
from galEixy formation physics and (in the case of this study) 
fi'om photo- 2; errors. 

For satellite galaxies, a power-law in mass provides a 
good description for the mean occupation number in simu- 
lations: 



a "cut-o 
V Mo j 



(27) 



Introducing a "cut-off" to the satellite HOD of the form 

'M-Mi\'^ 



< NalM >-- 



(28) 



(Zheng et al. 2005) did not significantly improve the fit of 
the model to the data. Our halo model power spectrum is 
hence specified by four variables: Mcut, fcut. Mo and /3. 

In addition we must satisfy one extra constraint: the 
galaxy rmmber density 
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dMn{M)N{M) 



(29) 



must match the observed rmmber density in each redshift 
slice. We match this constraint by fixing the variable Mcut 
for each choice of CTcut, Mo and /3. Our model hence contains 
three independent parameters. Figure 5 displays an example 
halo model correlation function and the component 1-halo 
and 2-halo terms. 

Useful quantities which can be derived from the HOD 
are the effective large-scale bias 



bg = I dMn{M)b{M) 



N{M) 



(30) 
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Figure 5. An example halo model galaxy correlation function at 
z = 0.5 for parameters Ucut = 0.3, logio(Mo/?i-iM-0) = 14.3 
and /? = 1.5. Matching galaxy density rig = 5 X 10~^ Mpc~^ 
requires logiQ{Mcut /h^^ Mq) = 12.98. The separate 1-halo and 
2-halo components arc shown, which cross-over in dominance at 
scale ~ lh~^ Mpc. This model has a linear bias bg = 1.82 and an 
effective halo mass \og'^Q{Mgg/h~^MQ) = 13.48. 



and the effective mass M^g of the halo occupation distribu- 
tion: 



Meff = / dM M n{M) 



N{M) 



(31) 



We can also determine the average fraction of central or 
satellite galaxies in the sample, e.g. for central galaxies: 



fc = 



J dMn{M)Nc{M) 



(32) 



/ dMn{M) Nc{M) [1 -I- NsiM)] 
and fs = i — fc for satellite galaxies. 

4.6 Conversion to angular correlation function 

Knowing the redshift distribution p{z) of a galaxy popula- 
tion, wc can project the spatial galaxy correlation function 
^(r) to an angular correlation function w{d) using Limber's 
equation: 



w{e) 



2 1 dxfixf du^{r= a/m2 _^a;26)2) (33) 
/o Jo 



where f{x) describes the radial distribution of sources as 
p{z) 



(34) 



dx/dz{z) 

where z is the redshift corresponding to co-moving radial 
co-ordinate x{z). 

5 PARAMETER FITS 

We fitted the 3-parameter halo model (fcut. Mo, /3) to the 
observed angular correlation functions in each redshift slice, 
fixing the remaining parameter Mcut by matching to the ob- 
served galaxy number density rig . We used a combination of 
a coarse grid-based search and a downhill-simplex method 
to locate the minimum value of the statistic, using the 
full covariance matrix. We then employed a fine grid-based 
method to explore the surface around the minimum and 



determine the errors in the fitted parameters (by marginal- 
izing over the other model parameters) . The mean and stan- 
dard deviation of each model parameter, marginalizing over 
the other parameters, are listed in Table 2. 

In Figure 6 we plot the best-fitting halo model correla- 
tion functions together with the data. We divide the results 
by the best-fitting power-law model from Section 3.3 for in- 
creased clarity. We note immediately that the halo model 
framework has successfully reproduced the deviations from 
the power-law, owing to the separate contributions of the 
1-halo and 2-halo terms. The result is a good fit of model 
to data, with the addition of only one extra parameter com- 
pared to the original power-law fit. The minimum values of 
are around 30 (for 24 degrees of freedom). Inspection 
of Figure 6 reveals that the best-fitting model always de- 
scribes the transition region between the 1-halo and 2-halo 
terms very well, and that the main source of discrepancy is 
at very small scales (< 200 kpc) where the data points lie 
systematically above the model prediction. 

In Table 2 we also list the derived values of the galaxy 
bias factor bg and effective halo mass Meff for each redshift 
slice (calculated using equations 30 and 31). The errors in 
these quantities were obtained by evaluating their values 
over the fine grid of halo model parameters and weighting 
by the appropriate probability for the fit of model to data 
at each grid point. The range of linear bias, increasing with 
redshift from bg ^ 1.92 (at z = 0.475) to bg ^ 2.62 (at 
z = 0.625), agrees well with fits to the large-scale power 
spectrum of the LRGs (Blake et al. 2007), allowing for the 
diflFering values of the normalization as (in this study we 
have assumed as = 0.8, whereas the bias values in Blake et 
al. (2007) are quoted for ag — 1). The effective mass, rang- 
ing from Meff = ioi»-6i^i3.8() /j-i j^j^ gygj. ^jjg gg^jj^g redshift 

range, confirms that LRGs are hosted by massive dark mat- 
ter haloes and are highly biased tracers of the clustering 
pattern. Our values for the effective mass are in the same 
range as the LRG halo mass measured by Mandelbaum et al. 
(2006b) using galaxy-galaxy lensing (M = iQ^-i ss j^-i 
for the bright sample of Mandelbaum et al.). We note that 
the systematic increase in galaxy bias in each redshift slice 
is driven more by the increasing luminosity threshold rather 
than redshift evolution. 

Figure 7 plots the statistical range of allowed halo oc- 
cupation distributions < A'^IM > in each of the four redshift 
slices assuming the parametric description of equations 26 
and 27. The parameters CTcut, Mo and j3 were varied over a 
grid and the probability determined at each grid point us- 
ing the x^ statistic. This probability distribution was used 
to construct the mean and standard deviation of the value 
of < A'^IM > as a function of halo mass M. The aver- 
age number of our galaxy sample hosted by a halo of mass 
M = 10^* '^ h-^ Mq is (5.5, 4.1, 2.6, 1.4) in the four redshift 
slices, decreasing systematically with redshift as the thresh- 
old luminosity increases. Broadly speaking, the effect of in- 
creasing luminosity is to shift the HOD to higher masses 
without significantly changing its shape, i.e. shifting to the 
right in Figure 7 in a similar fashion for central and satel- 
lite galaxies. If we weight the HOD by the mass function 
of haloes which steeply decreases with increasing mass, we 
find that the fraction of galaxies that are classified as cen- 
tral galaxies is very high: /c = (0.88, 0.90, 0.93, 0.97) in the 
four redshift slices. This agrees with the standard picture 
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Figure 6. The measured angular correlation functions in the four redshift slices together with the best-fitting halo models. The y-axis is 

the ratio of the angular correlation function w{6) and the best-fitting power-law model in each slice (from Section 3.3). The co-moving 
spatial scale corresponding to angular separation 0.1° at the median redshift of each slice is displayed in the bottom left-hand corner of 
each plot. A flat ACDM Universe is assumed with cosmological parameters Qm = 0.3, h = 0.7, Hb/flni = 0.15, as = 0.8 and Ub = 1- 



Table 2. Halo model fits to the angular correlation function in the four redshift slices. Column 2 records the observed galaxy number 
density which our models are constrained to match. Columns 3, 4 and 5 list the best-fitting values and errors of the halo model parameters 
Ccut , Afo and f3 defined by equations 26 and 27. Column 6 displays the inferred value of the final halo model parameter Mcut . In Columns 
7 and 8 we evaluate the corresponding galaxy bias factor bg (equation 30) and effective halo mass M^g (equation 31). Column 9 records 
the minimum value of the chi-squaxed statistic for the halo model, Xj^aio' evaluated using the full covariance matrix (which we compare 
with 24 degrees of freedom). Base- 10 logarithms are used in this Table. 



Redshift slice 10^ x n, acut log (l^) (l^) (i^) ^h, 

(h^ Mpc-3) 



0.45 < z <0.5 


5.03 


0.21 ± 0.11 


14.09 ±0.01 


1.57 ±0.02 


12.98 


1.92 ±0.02 


13.61 ±0.01 


22.1 


0.5 < z < 0.55 


3.07 


0.07 ±0.07 


14.22 ±0.01 


1.69 ±0.04 


13.12 


2.15 ±0.01 


13.67 ±0.01 


28.3 


0.55 < 2 < 0.6 


1.60 


0.24 ±0.12 


14.39 ±0.01 


1.87 ±0.07 


13.35 


2.38 ±0.04 


13.74 ±0.01 


36.8 


0.6 <z < 0.65 


0.56 


0.53 ±0.14 


14.76 ±0.10 


1.80 ±0.36 


13.79 


2.62 ± 0.09 


13.80 ± 0.03 


40.1 



of Luminous Red Galcixies forming at the heart of massive 
dark matter haloes. 

The best-fitting value of /3, the slope of the power-law 
HOD for satellite galaxies, increases systematically with red- 
shift. The most important cause of this trend is evolution 
of /3 with luminosity, as measured at low redshift (Zehavi et 
al. 2005b). Investigating this further, we re-fitted the halo 
model parameters in each of the redshift slices for each of 
the absolute i-band magnitude thresholds Mi — Slog^g/i = 
(-22.23, -22.56, -22.87, -23.20) which correspond to the 
luminosity thresholds of the four redshift slices. Hence all 



redshift slices can provide a sample with Mi — blogi^h < 
—23.20, but only the first redshift slice contributes a sample 
with Mi - 51ogio/i < -22.23. We obtained the probabil- 
ity distribution for /? for each halo model fit, marginalizing 
over the other model parameters, and combined the results 
for matched luminosity samples in different redshift slices. 
We find that the measurements of /? at a fixed luminosity 
threshold agree well between the redshift slices, and the com- 
bined result for the four luminosity thresholds listed above 
is 13 = (1.57±0.02, 1.68±0.03, 1.91 ±0.05, 2.15 ±0.14), con- 
firming a significant evolution of /? with luminosity. Consid- 
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Figure 7. The range of halo occupation distributions fitted to 
each of the four redshift slices (with each subsequent slice shifting 
from left to right in the Figure) assuming the parametric descrip- 
tion of equations 26 and 27. This Figure was generated by varying 
the parameters (Tcut, Mq and /3 over a grid and determining the 
relative probability at each grid point using the statistic. At 
each halo mass M, the probability distribution within the model 
parameter space was used to construct the mean HOD (the solid 
line) and the 68% confidence region (the dotted lines). 

ering just the first redshift slice, the best-fitting effective halo 
mass for the four luminosity slices is \ogiQ{Mcg /h~^ Mq) = 
(13.61, 13.73, 13.88, 14.02) corresponding to number densi- 
ties Ug = (5.03, 2.45, 1.02, 0.32) x 10"* Mpc~^ 

Our results have a significant dependence on the as- 
sumed value of (Tg = 0.8, which sets the overall normal- 
ization of the matter power spectrum. For erg — 0.7, the 
best-fitting effective halo mass in each of the four redshift 
slices is logio(Mcff//i-^Mo) = (13.48,13.54,13.61,13.68). 
For erg = 0.9, the results were (13.72,13.79,13.85,13.91). 
In all cases the minimum values of were similar to those 
obtained for our default value of erg, suggesting that there 
is a strong degeneracy between erg and the halo model pa- 
rameters used in this analysis. The changing ag affected the 
best-fitting value of Mo much more strongly than the value 
of /3. We also tried lowering the value of the scalar spec- 
tral index of the primordial power spectrum from Ws = 1 to 
Us = 0.95, as supported by recent observations of the Cos- 
mic Microwave Background radiation (Spergel et al. 2007). 
Assuming as = 0.8, the best-fitting effective masses are 
(13.59, 13.66, 13.72, 13.78) which do not differ significantly 
from our default values presented in Table 2. 



6 COMPARISON TO PREVIOUS WORK 

Several previous studies have fitted halo model parameters 
to populations of red galaxies, for example: Magliocchetti & 
Porciani (2003, 2dFGRS); Zehavi et al. (2004, 2005b, SDSS); 
Collister & Lahav (2005, 2dFGRS groups catalogue); Phleps 
et al. (2006, COMBO-17 survey); and White et al. (2007, 
NDWFS). We make comparisons to these analyses below, 
where possible. 

Halo model fits to the 2dFGRS galaxy correlation func- 
tion for late-type and early-type galaxies were performed by 
Magliocchetti & Porciani (2003). In addition, Collister & La- 
hav (2005) directly investigated the distribution of galaxies 



within 2dFGRS groups. These two studies produced a rea- 
sonably consistent measurement of the slope /? ~ 1 of the 
HOD at high masses for early-type galaxies. Our best-fitting 
slope, /3 = 1.5 2.0, is much higher due to two factors: (1) 
the significantly higher luminosity of our galaxy samples; (2) 
we make the distinction between central and satellite galax- 
ies, separating out a central galaxy contribution A'^c ~ 1 at 
high masses. This latter has the effect of significantly steep- 
ening the slope of the power-law HOD fitted to the remain- 
ing satellite galaxies (which in fact contribute only 5 — 10% 
of our sample, as noted in Section 5). In other words, we 
effectively fit a model A^ = Nc{l + A''^) ^ 1 + {M/Mof at 
high masses, rather than N = (M/Mo)''. 

Zehavi et al. (2004) analyzed a luminous subset of 
galaxies from the SDSS "main" spectroscopic database with 
Mr < —21 and mean redshift z ~ 0.1. They found that a 
HOD of the form 

< N\M > = (M < Mcut) 

= 1 (Mcut < M < Mo) 

= (M/Mo)" {M > Mo) (35) 

produced a good fit to the clustering data, where Mcut — 
^q12.79^-i^^^ Mo = lO"-^8ft-^M0 and /3 = 0.89. The 
effective mass corresponding to these parameters is Mes = 
lO" *^/i"^M0 (for their choice of erg = 0.9). The number 
density of the Zehavi et al. (2004) sample is rig = 9.9 x 
10~*/i^ Mpc~^ which is a factor of 2 — 3 higher than our 
study. Although the effective masses are similar, we note 
that the redshift difference between the Zehavi et al. (2004) 
sample and ours may be important; the number density of 
a sample of dark matter haloes of fixed mass increases with 
decreasing redshift owing to the growth of structure. The 
difference in the best-fitting value of the power-law slope 
between Zehavi et al. and our analysis is again connected 
to the different forms of HOD fitted (owing to our central 
galaxy contribution, as discussed above). 

Zehavi et al. (2005b) presented an extended analysis of 
the SDSS data in which central and satellite galaxy con- 
tributions are considered separately. Their default model 
includes a sharp cut-off for the central galaxy HOD at a 
fixed mass, rather than our "softened" transition from to 
1 galaxies. They find that the slope of the power-law satellite 
HOD increases systematically with luminosity in a manner 
entirely consistent with our high- luminosity measurements 
of /? = 1.5 ^ 2.0. In addition, Zehavi et al. (2005b) note 
that the step function for < A^dM > produces a poor fit 
to the data in their highest-luminosity bin, consistent with 
our requirement for a softened transition parameterized by 
(Jcut. They also find, in agreement with our analysis, that 
the great majority of luminous galaxies are central galaxies 
of their host dark matter haloes, rather than satellites in 
more massive systems. A low ( < 10%) satellite fraction for 
the most luminous elliptical galaxies is also found in galaxy- 
galaxy lensing studies (Seljak et al. 2005; Mandelbaum et 
al. 2006a) and other clustering studies (Tinker et al. 2007; 
van den Bosch et al. 2007). 

Phleps et al. studied various populations of galaxies in 
the COMBO-17 survey at a mean redshift J = 0.6 which is 
similar to our dataset. For red-sequence galaxies, Phleps et 
al. quote an effective halo mass for their best-fitting model of 
A'l^g = lO^^'^ h-^MQ, whereas we find Afcff = 10^^-'' K'^Mq 
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(Table 2). This apparently large discrepancy is caused by the 
significant difference in the luminosity threshold of the two 
samples: the number density of our LRG catalogue is more 
than an order of magnitude smaller (there is also a difference 
in the assumed value of ag). 

White et al. (2007) fitted a Halo Occupation Distri- 
bution model to the clustering of Luminous Red Galaxies 
in the NOAO Deep Wide-Field Survey (NDWFS) Bootes 
field of 9 dog'^, analyzed in redshift slices between z = 0.4 
and z = 1.0. The luminosity thresholds arc fixed such that 
the galaxy number density in each redshift slice is 10~^ 
Mpc~^, exceeding our sample by a factor w 3 at a = 0.5. 
White ct al. demonstrated that the clustering of the z — 0.5 
sample cannot be accounted for by simple passive evolution 
of the z = 0.9 sample, but rather there must be merging 
or disruption of the most luminous satellite galaxies in mas- 
sive haloes. The best-fitting satellite fraction in the NDWFS 
sample is found to be 18%, a little higher than the results 
of our study, but consistent with a trend in which satellite 
fraction decreases with increasing luminosity. 

In conclusion, our halo model parameter measurements 
appear broadly consistent with previous work, allowing for 
differing luminosity thresholds. A fully consistent compari- 
son of our analysis at z ~ 0.55 with results at 2 ~ is beyond 
the scope of this work, owing to the differing forms of halo 
occupation distribution assumed by different authors, but a 
topic worthy of further investigation. 

Measurement of the 3-point clustering functions will 
add further insight into the LRG clustering properties. Re- 
cent work by Kulkaxni et aJ. (2007), analyzing the SDSS 
spectroscopic LRG sample at 2 w 0.35, favoured a shallower 
slope for the satellite HOD, j3 w 1.4, with a higher satel- 
lite fraction of 17%. Further study is required to understand 
these differences. 



7 TESTS FOR SYSTEMATIC PHOTOMETRIC 
ERRORS 

We performed a series of tests for potential systematic pho- 
tometric errors that may affect our clustering results. Fol- 
lowing the discussion in Blake et al. (2007), we compared 
the angular correlation function measured for the "default" 
sample with that obtained by restricting or extending the 
galaxy selection in the following ways: 

• Exclusion of areas of high dust extinction (> 0.1 mag). 

• Exclusion of areas of poor astronomical seeing (> 1.5 

arcscc) . 

• Exclusion of areas lying in the overlap regions between 
survey stripes. 

• Exclusion of areas in the vicinity of very bright objects 
(circular masks of radius 1 arcmin around objects with i < 
12). 

• Variations in the star-galaxy separation criteria. This 
is quantified by the coefficient Sag in Blake et al. (2007) and 
CoUister et al. (2007), which fixes the aggressiveness of the 
star-galaxy separation in the neural network. Our default 
choice is S^g > 0.2, which results in a level of stellar contam- 
ination of 1.5% with the loss of only ~ 0.1% of the genuine 
galaxies (see Fig. 13 in Collister et al. 2007). We also tried 
Ssg > (no additional star-galaxy separation in the neural 



network; stellar contamination 4.4%) and 5sg > 0.8 (stellar 
contamination 0.4%; loss of 1.2% genuine galaxies). 

We refer the reader to Blake et al. (2007) for a more thor- 
ough discussion of the possible effects of these systematic 
errors. 

Our results are presented in Figures 8, 9 and 10. Each 
plot is composed of four panels, one for each redshift slice. In 
each panel we show the angular correlation function for the 
default sample together with that corresponding to a change 
in the galaxy selection criteria. We divide all the correlation 
functions by a power-law fit to the default measurement to 
render the results more easily comparable. 

We conclude from Figures 8 and 9 that our results are 
robust against the details of the angular selection function: 
varying dust extinction, seeing, overlap regions and bright 
object masks all have little effect on the measured corre- 
lation function. Figure 10 reveals that the details of the 
star-galaxy separation affect the amplitude of the measured 
correlation function although not (to first order) the shape. 
This amplitude shift is already encoded in the stellar con- 
tamination factor (1 — Z)'^. In no case does a change in the 
galaxy selection alter the detectability or shape of the halo 
model signature. 



8 CONCLUSIONS 

We have measured the angular correlation function of LRGs 

in the SDSS imaging survey, using accurate photometric red- 
shifts to divide the galaxies into narrow redshift slices and 
create volume-limited samples. We find that: 

• A canonical power-law fit provides a poor description 
of the small-scale angular correlation function, although the 
best-fitting slope w{9) oc 0~°-^^ agrees well with previous 
studies of Luminous Red Galaxies. 

• The halo model of galaxy clustering, composed of sepa- 
rate 1-haIo and 2-halo contributions, produces a good fit to 
the deviations from a power-law. We assume a halo occupa- 
tion distribution with separate components for central and 
satellite galaxies, implementing realistic models for scale- 
dependent bias, halo exclusion and non-linear growth of 
structure. We find that the HOD for central galaxies re- 
quires a "soft" transition from zero to one central galaxy, as 
opposed to a step function, to re-produce the observations. 
The functional form A''^ = 0.5{1 -|- erf [logio(M/Mcut)/o-cut]} 
provides a good fit for central galaxies, combined with a 
power-law HOD for satellite galaxies = (M/Mo)". One 
parameter of the model (Afcut) is fixed by the overall number 
density of the galaxy sample; hence this halo model contains 
3 variable parameters (acut. Mo and /3), just one more than 
a simple power law. Allowing the concentration parameter c 
to vary, or including a more sophisticated HOD for satellite 
galaxies, does not improve the model fits. 

• The slope /3 of the power-law HOD for satellite galaxies 
is a strong function of luminosity, increasing to /3 « 2 for our 
most luminous sample. This is consistent with extrapolating 
the variation of /3 with luminosity in local samples (Zehavi 
ct al. 2005b). The physical implication of this result is that 
haloes of higher mass have greater relative efficiency at pro- 
ducing high-luminosity satellites. We find no variation of f3 
with redshift across our sample (from z = 0.45 to z = 0.65). 
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Variations in angular selection function (1) 



0.45 < z < 0.50 



Default 

Extinction < 0.1 mag 
Seeing < 1 .5 orcsec 




All data divided by povver-law fit to Defoult 
0,01 0.1 



O.dO < z < 0.55 



Default 

Extinction < 0.1 mag 
Seeing < 1 .5 orcsec 



1 =1o- 




• Default 

o Extinction < 0.1 mag 
A Seeing < 1 .5 arcsec 




/ deg 

Figure 8. The dependence of the angular correlation function measurement in 4 redshift slices on varying dust extinction and astronomical 
seeing. Results from the default catalogue are compared to an analysis restricting the regions analyzed to (1) a maximum dust extinction 
of 0.1 mag, or (2) a maximum seeing of 1.5 arcsec. The j/-axis is the angular correlation function w{0) divided by a power-law fit to the 
default model. 



Variations in angular selection function (2) 



0.45 < z < 0.50 



Default 

Star masks < 12 mag 
Excluding overlop regions 




All data divided by power-law fit to Default 
0.01 0.1 



O.dO < z < 0.55 



Default 

Star miQsks < 12 mag 
Excluding overlop regions 




0,;o < z < 0,60 



<I1 



• Default 

o Star masks < 12 mag 
A Excluding overlap regions 




• Default 

o Star masks < 12 mag 
ii Excluding overlap regions 




B / deg 

Figure 9. The dependence of the angular correlation function measurement in 4 redshift slices on the presence of bright objects and 
stripe overlap regions. Results from the default catalogue are compared to (1) an analysis placing circular masks of radius 1 arcmin 
around all objects with i-band magnitudes brighter than 12, and (2) an analysis excluding overlap regions between stripes. The results 
are displayed in the same manner as Figure 8. 
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Variations in star— galaxy separation 



0.45 < z < 0.50 



• Default 

More strict {6^^ = 0.8) 

1 More liberal {6 = 0.0) 



All data divided by power-law fit to Defoult 



2 "^o"' 



1 So- 



.br> < z < 0.60 



• Default 

o More strict (S^^ = 0.8) 
A More liberal (5 = 0.0) 



• Default 

o More strict (6^^ = 0.8) 
A More liberal (6 = 0.0) 




B / deg 

Figure 10. The dependence of the angular correlation function measurement in 4 redshift slices on the star-galaxy separation criteria. 
Results from the default catalogue are compared to analyses changing the value of the co-efRcient 5sg defined in Blake et al. (2007) and 
CoUister et al. (2007), which controls the aggressiveness of the star-galaxy separation. The results are displayed in the same manner as 
Figure 8. 



• The best-fitting width a^ut of the transition from zero 
to one central galaxy with increasing mass is in the range 
CTcut = 0.1 — > 0.5 (this is the parameter least- well con- 
strained by our analysis). The physical implication of this 
result is a scatter between central galaxy luminosity and 
host halo mass. This scatter results from galaxy formation 
physics and from the photo-z error in the conversion of ap- 
parent magnitude to luminosity. 

• The halo model fits describe how Luminous Red Galax- 
ies populate dark matter haloes as a function of their mass 
M. The average number of our galaxy sample hosted by a 
halo of mass M = lO"'^ /i"^ Mq is (5.5,4.1, 2.6, 1.4) in the 
four redshift slices, decreasing systematically with redshift 
as the threshold luminosity increases. Broadly speaking, the 
effect of increasing luminosity is to shift the HOD uniformly 
to higher masses without significantly changing its shape. 
The large majority of galaxies in our sample are classified 
as central galaxies of their host dark matter haloes, rather 
than satellites in more massive systems, in agreement with 
previous studies of galaxy-galaxy lensing and clustering of 
the most luminous galaxies. The satellite fraction varies in 
the range 3% — »■ 12% across the redshift slices. 

• The halo model fits provide robust predictions of the 
average linear bias of the LRGs on large scales and the 
effective mass of their host dark matter haloes. The re- 
sulting amplitude of the linear bias (6g = 1.9 — > 2.6, 
increasing with redshift, assuming a normalization of the 
matter power spectrum as = 0.8) agrees well with fits to 
the large-scale power spectrum (Blake et al. 2007). The 
effective halo mass {M^a = lO^^'^^^^-* /i"^Mo) provides 
a quantitative statement of how the LRGs trace the un- 



derlying dark matter haloes. The value of Mes has a sig- 
nificant dependence on as, and a weak dependence on 
the slope of the primordial scalar index rig: the effective 
mass increases by Alog^g(Meff //i~^Mq) w 0.1 when the 
value of (78 is increased from 0.8 to 0.9, and decreases by 
Alogio(Mefi//i"^M0) « 0.02 when Us is decreased from 1.0 
to 0.95. 

Future studies will explore joint fits of the cosmological 
parameters and halo model parameters (Abazajian et al. 
2005; Zheng & Weinberg 2007); direct measurement of halo 

occupation via a cluster and group-finding analysis of the 
photo-z catalogue; a consistent comparison of clustering of 
Luminous Red Galaxies at z « 0.5 and at z « 0; and testing 
the halo model further via 3-point clustering statistics and 
higher moments. 
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